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energy range relevant for contained neutrino interactions, and we identify the major 



We compare several different calculations of the atmospheric neutrino flux in the 




sources of difference among the calculations. We find nothing that would affect the 



predicted ratio of Vel^pu^ which is nearly the same in all calculations. Significant dif- 



ferences in normalization arise primarily from different treatment of pion production 



> 




by interactions of protons in the atmosphere. Different assumptions about the pri- 



mary spectrum and treatment of the geomagnetic field are also of some importance. 
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I. INTRODUCTION 



Two deep underground detectors [|T],3 observe a ratio of electrons to muons that is signifi- 
cantly different from what is expected from the spectra of z/g and z/^ in the atmosphere. These 
two experiments use large volumes of water to detect Cherenkov radiation from charged par- 
ticles that originate from interactions of neutrinos inside the detector. Together these two 
experiments have collected about 80% of the world's statistics of atmospheric neutrino inter- 
actions. Events with a single Cherenkov ring, which are mostly quasi-elastic, charged-current 
interactions of neutrinos, constitute the simplest and largest class of events in these detec- 
tors. The observed ratio of electron-like to muon-like events is significantly greater than 
expected from calculations. 

A major component of the calculations (along with neutrino cross-sections and detector 
response) is the evaluation of fluxes of atmospheric neutrinos. Four sets of atmospheric 
neutrino spectra have been used in the past several years to interpret the measurements of 
interactions of z/g (z/g) and z/^ (z/^) in underground detectors. All four flux calculations agree 
within a range of 5% for the flavor ratio of neutrinos with 0.4 < < 1 GeV 0, which is 
perhaps not surprising as most of the sources of uncertainty cancel in the calculation of this 
ratio. Much larger differences exist among the results for normalization and shape of the 
spectra, and these lead to ambiguities in the interpretation of the anomalous flavor ratio. For 
example, the calculation of Bugaev and Naumov (BN) [Q has a harder spectrum than the 
other calculations (relatively less neutrinos below E^, = 500 MeV than above 1 GeV). 

Such a hard spectrum allows the suggestion |^ that the observed relative excess of electron- 
like events could be the result of the decay p e^vv. Even if spectra have the same shape, 
overall differences in normalization suggest different interpretations of the anomaly in terms 
of neutrino oscillations. A high normalization that agrees with the observed electron flux 
favors oscillations predominantly in the z/^ ^ v^- sector whereas a low normalization 
would suggest an oscillation that includes z/g. A calculation with a low normalization relative 
to which there is an excess of electrons could also be consistent with an interpretation in 
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terms of a neutron-induced background masquerading as interactions of z/g [|T0| . 



Of the four calculations we consider, three are completely independent of each other. 
The independent calculations are Honda et al. (HKHM) [Q, Bugaev & Naumov (BN) |^ 
and Barr, Gaisser & Stanev (BGS) 0. The work by Lee & Koh (LK) ^ uses a three- 
dimensional version of the model of hadronic interactions from BGS. LK also use the same 
primary spectrum as BGS. We have discovered several bugs in the implementation of the 
LK code. When these are removed, the results of LK are essentially the same as those of 
BGS in the absence of a geomagnetic cutoff. Moreover, since the transverse momentum of 
a neutrino from decay of a pion or muon is typically no more than 30 MeV, for the energies 



of interest here {E^, > 200 MeV) a three dimensional calculation is not necessary |TT|,[T2|. In 
what follows we therefore do not consider separately the calculation of LK. 

The calculation of BGS is a Monte Carlo simulation made in two steps: first, cascades 
are generated for primary protons at a series of fixed primary energies over an appropriate 
range of zenith angles; second, the resulting yields of neutrinos, binned in E^, are added 
together after weighting by the primary spectrum and geomagnetic cutoffs characteristic of 



each detector location, as calculated in Ref . |]TE[ , which neglected the effect of the penumbra 



of the Earth. This structure allows us to substitute other assumptions about primary spectra 
and composition and about geomagnetic cutoffs while keeping all other input the unchanged. 
Thus we can compare the sensitivity to these assumptions one-by-one in isolation. Hadron 
production in BGS is described in a model called TARGET [|I4[ which is a parametrization of 



accelerator data for hadron nucleus collisions with emphasis on interaction energies around 
20 GeV, which are most important for production of GeV neutrinos. 

The calculation of BN is a semi-analytic integration of the atmospheric cascade equations 
in "straight-forward approximation" over the primary spectrum as modified by appropriate 
geomagnetic cutoff rigidities for protons and nuclei. Namely, BN used detailed tables for 
effective vertical cutoffs from Ref. |jl5[ (corrected for the displacement of the geomagnetic 
poles with time) and a dipole-like relation for the directions different from vertical. In this 
approach, the penumbra structure, contribution of re-entrant albedo and direct influence of 



the geomagnetic field on the charged secondaries were neg lectedS. It has been estimated ITOll 
that these effects are at the level of 10% or less for atmospheric neutrino fluxes averaged 
over reasonably wide solid angle bins. 

The hadronic interaction model used by BN is an analytic parametrization of double- 
differential inclusive cross sections [ITS], which is based on great array of accelerator data 



and, according to [|I^, it applies at > 1 GeV/c, > 150 MeV/c and pk± > 300 
MeV/c. The comparison of the model with some data, which was not used when fitting its 
parameters, was presented in Ref. [|l6l (see also Ref. [l^])- The exact inclusive kinematics 
was drawn on to make all necessary integrations. Other details of the BN calculation were 



described in Refs. [19,20,16 



The work of HKHM is a Monte Carlo calculation that includes a detailed treatment 
of the effect of the geomagnetic field. Cutoffs are calculated for each detector location by 
backtracing antiprotons through a map of the geomagnetic field. This procedure was also 



used by LK, and a similar analysis has recently been carried out by Lipari and Stanev [21 



This is the correct way to account for cutoffs because it includes the effects of trajectories 
that are forbidden because they intersect the surface of the Earth (penumbra). For the 
interactions above 5 GeV, HKHM use the subroutine packages FRITIOF version 1.6 [122 
and JETSET version 6.3 H. At lower energy the algorithm NUCRIN [|2^ is used. 

All calculations include the effect of muon polarization on the neutrinos from muon decay, 
following the remark of Volkova [^] who emphasized its importance in this context!. 

A quantitative comparison of the three independent calculations is made in Table 1. 
The first two blocks show the neutrino fluxes (normalized to BGS) in three ranges of energy, 



^The last two effects were also neglected in the other calculations. 

^The fluxes shown by BN in Ref. [^ do not include muon polarization, but they have since been 
corrected for this effect together with the effect of muon depolarization caused by muon energy 



loss [26, l"; 
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0.4 < E,y < 1, 1 < Ely < 2, and 2 < < ?> GeV. The third block compares the neutrino 
ratios in the energy range 0.4 < E^ < \ GeV. We tabulate Re/^i = {i^e + li^e)/i^fj. + ^i^/j.) 
because the cross section for quasi-elastic interactions of antineutrinos is roughly one-third 
that of neutrinos in the low energy range relevant for single-ring contained events. Table 1 
demonstrates the main differences between the different calculations. In addition to the 
overall normalization, there are also significant differences in the z/g/ z/g ratio. The difference 
between BGS and HKHM, 0.89 vs. 0.84 averaged over the 0.4 to 1 GeV range, is actually 
quite big at neutrino energy below 500 MeV and disappears above 1 GeV. 

We divide our discussion into three sections. We first compare the assumptions about 
the primary spectrum and about the geomagnetic cutoffs made in the three calculations. We 
then compare the treatment of hadronic interactions. We conclude with a brief discussion of 
the implications for interpretation of the measurements of contained neutrino interactions 
and the anomalous flavor ratio of neutrinos. 

II. PRIMARY SPECTRUM, COMPOSITION AND GEOMAGNETIC CUTOFFS 

The cosmic ray beam consists of free protons and of nucleons bound in nuclei. The two 
must be distinguished because the primary spectrum and geomagnetic cutoffs depend on 
rigidity, while production of secondaries depends on energy per nucleon. At fixed rigidity, 
the momentum per nucleon for nuclei is half that for free protons to a good approximation. 

BGS assume equal numbers of neutrons and protons for the bound nucleons and obtain 
the neutron yields from the proton yields by reversing all electric charges of mesons. This 
is exact for pions (by isospin symmetry) but introduces a small excess of charged kaons 
in interactions of neutrons. However, kaons contribute less than 10% of the neutrinos in 
the energy range relevant for contained events (E. Frank, private communication, 1993). A 
separate calculation for incident protons and neutrons has since been made |12| that has 
verified that this approximation has a negligible effect for the energy range relevant for 
contained interactions. 
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All bound nucleons are assumed to interact with the same starting point distribution 
as free nucleons; i.e., their cascades are simulated exactly as if they were free. This as- 
sumption can be derived for neutrinos (where correlations between different nucleons in the 
nucleus are irrelevant) within the framework of the Glauber multiple scattering picture of 
nucleus-nucleus collisions [P7|. In contrast, BN treat nuclear interactions separately using 
the nucleus-nucleus cross sections to determine where the interactions occur |T9|j2( 

Given the two-part structure of the BGS calculation, it is possible to check the effect of 
the different treatment of primary spectrum (including the model of nuclear interactions) 
separately from all other effects. We have done this by using the BN and the HKHM 
primary spectra and the BGS yields and geomagnetic cutoffs. Table 2 shows the effective 
differential spectrum of primary nucleons for BGS and BN at solar minimum. Replacing the 
BGS primary spectrum with that of BN decreases the neutrino flux by only about 5% below 
1 GeV, leaves it nearly unchanged between 1 and 2 GeV and increases the calculated flux 
by about 5% for 2 < ii^j, < 3 GeV. The neutrino flavor ratio also remains largely unchanged. 
Thus the differences in the primary spectrum and the treatment of nuclear projectiles are 
not the origin of the significant differences between these two calculations. 

The primary spectra of BGS and HKHM averaged over the solar cycle are compared in 
Table 3. Use of the HKHM primary spectrum instead of that of BGS increases the neutrino 
flux in the energy range 0.4 - 1 GeV by about 7% and by about 12% in the 1-2 GeV 
range. The flux ratios again remain unchanged. In addition the BGS primary spectrum is 
significantly heavier, i.e. contains more bound nucleons, than either BN or HKHM spectra. 

To check for the influence of the geomagnetic cutoff models, we replace the offset dipole 
model (with no penumbra) [|13| used by BGS with the detailed treatment of the cutoffs used 
by HKHM. The flux of neutrinos between 0.4 and 1 GeV decreases by 10%. The geomagnetic 
effect is less important at higher energy. Note that the direction of this change is opposite to 
that caused by the different primary spectra used by BGS and HKHM, which was discussed 
above. We discuss this point further in the concluding section. 
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III. TREATMENT OF HADRONIC INTERACTIONS 



A useful way to compare the interaction models is to evaluate the spectrum-weighted 
moments (Z-factors) of the inclusive cross sections in the energy range relevant for the 
contained neutrino events. The Z-factors were not actually used in any of the calculations 
but they characterize the efficiency of the particle physics model to produce secondary 
particles in atmospheric cascades. The most important range of interaction energies for 
production of neutrinos with energies from 300 MeV to 3 GeV is ~ 5 < En < 50 GeV ||30|| . 
In this energy range cross sections do not scale so the Z-factors are energy-dependent: 



where x = Et,/Ej^, Ej^ is the total energy of the incident nucleon in the lab system and E^^ 
is the energy of the produced pion. 

The Z-factors are shown in Figure 1. The yields are rather different in the three sets 
of models. The decrease of Z^ch at incident energies above 10 GeV in the calculations of 
HKHM and BGS is correlated with the onset of kaon production. The curves of the 
BN calculation |^ have a different behaviour, reflecting less production of low energy pions 
(see Fig. 5 below). This key difference is directly related to the hardness of the BN neutrino 
energy spectrum in Table 1. 

In addition to the differences in energy- dependence and magnitude of the Z-factors, there 
are also some significant differences in the charge ratios. For almost all of the relevant range 
of energies, the ratio Zp^T^+ /Zp^^~ is significantly larger for BN than for either FRITIOF 
(HKHM En>5 GeV) or TARGET (BGS), as shown in Fig. 2. The 5 GeV value of NUCRIN 
(HKHM) is intermediate, and at lower energy NUCRIN gives a very high value of the ratio. 
This is an artifact of NUCRIN, which neglects the charge exchange process p + n +p, 
followed by A° p-K^. These differences are directly relevant to the differences in the ratio 
Re = Vej i^e the three calculations. 

To estimate the effect of these differences on the neutrino flux, we plot the yield of charged 
pions in three different momentum windows, weighted by the primary spectrum E^'"^ . This 
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is shown in Fig. 3 for pion momenta 0.3-0.4 GeV/c, 3-4 GeV/c and 6-8 GeV/c. The area 
under the curves is proportional to the pion flux in the three momentum windows, which 
are chosen to reflect three different characteristic neutrino energies: ~ 100 MeV, ~ 1 GeV 
and ~ 2 GeV respectively. The ratios BN/BGS for the three energy bands is respectively 
~ 0.3, ~ 0.6, and ~ 0.8 - 0.9. These roughly correspond to the ratios of BN/BGS in Table 
1. We therefore conclude that the difference in treatment of pion production in proton- air 
interactions is the main source of the difference between these two flux calculations. 

To explore further the differences between the interaction models we compare in Fig. 4 
the total multiplicity of charged pions in the three representations of hadronic interactions 



to data on proton-proton collisions [^. Both BGS and HKHM have significantly higher 
multiplicity than BN, which is similar to pp collisions. Both BGS and HKHM have a 
multiplicity on nuclear targets (A=14.5) that is about 50% higher than for proton targets. 
The pion production in the three interaction models is compared with data on light nuclei 
(beryllium) in Fig. 5. This data [0,^] is for beam momenta in the range 19-24 GeV/c, 
which is the median energy for production of ~ GeV neutrinos All three models fit 
the data well for x ~ 0.2. At smaller x BN has much lower pion yield than the other two 
models, which leads to a correspondingly low result for the calculated neutrino flux below 
1 GeV. The difference becomes progressively less with increasing neutrino energy because the 
representations of pion production agree with each other rather well at higher pion energy. 

Much of the available accelerator data in the relevant range of beam momentum for 
protons incident on light nuclei is from experiments carried out to calculate accelerator 



neutrino beams. For example, the Eichten et al. experiment was performed to provide 
data for calculation of the neutrino beam at the CERN-PS. A survey used at Brookhaven 
is that of Sanford & Wang p2l. Data used to calculate the Argonne neutrino beam is 



summarized by Barish et al. |3^. In all these cases the data are for relatively high pion 
momentum, so the ambiguity at low momenta in Fig. 5 is difficult to resolve. If we use data 



on proton targets at 24 GeV M^ to extend the data into the region of x < 0.15, then the 
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data would favor the model used by BN.0 We emphasize, however, that what is relevant is 
proton interactions in nitrogen and oxygen nuclei, in which the multiplicity of low energy 
pions is bound to be enhanced to some degree. The model of BGS and the Lund FRITIOF 
model as used by HKHM both give significantly higher pion yields at low x than that 
used by BN. We consider this question to be unresolved at present. 



IV. SUMMARY 

We find that differences in the representation of the production of pions in ~ 10 to 
~ 30 GeV interactions of protons with light nuclei are the major source of differences among 
three independent calculations of the flux of atmospheric neutrinos in the GeV energy range. 
Approximate treatment of the geomagnetic cutoffs also contributes significantly, especially 
for low energy neutrinos at Kamioka, which is the site with the highest geomagnetic cutoff 
for downward cosmic rays. 

The lower neutrino flux below 1 GeV in the BN calculation is due to the representation of 
pion production they use, which is similar to interactions on a nucleon target. On the other 
hand, both BGS and HKHM use representations of proton-nucleus data that correspond to 
a pion multiplicity enhanced by a factor 1.6 as compared to proton-proton collisions. 

The approximate representation of the geomagnetic cutoffs used by BGS tends to overes- 
timate the intensity of low energy cosmic rays that penetrate through the geomagnetic field 
to the atmosphere. This is mainly because the effect of the Earth's penumbra has been ne- 
glected. The preferred approach is to use the ray-tracing technique, as done by HKHM and 
LK. In a new calculation of the geomagnetic cutoffs, including all details properly, Lipari and 
Stanev |^ find neutrino fluxes at 1 GeV reduced by factors of 0.85, 0.92 and 0.97 relative to 



BGS respectively at Kamioka, Gran Sasso/Frejus and IMB/Soudan. An updated version of 



^We are grateful to D.H. Perkins for pointing out this reference to us and making this comparison 
51. 
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BGS neutrino fluxes , based on the correct treatment of tlie geomagnetic cutoffs, predicts 
neutrino fluxes equal to, or slightly lower than, those of HKHM. 

We have referred earlier to the fact that the primary spectrum used by HKHM is sig- 
nificantly higher than that of the other two calculations. In addition, it shows a very large 
variation between solar maximum and solar minimum, even at equatorial latitudes. On the 



other hand, the primary spectrum of BGS gives a solar cycle variation that is proba- 
bly too small. The level of uncertainties associated with the primary spectrum is currently 
±10%. The normalization of the primary cosmic ray spectrum and its dependence on the 
solar cycle epoch should be a subject of future studies. 

The differences in the Ve/i^e ratio, which are not yet a major factor in the interpretation 
of the experimental results, have to be further explored. 

Comparison to high altitude muons may be used to check the normalization in a global 
way that may avoid the need to resolve all the various differences in the input to a calculation 



that starts from the primary cosmic ray spectrum ||3^. There are recent measurements of 
muons at high altitude PB| to which calculations discussed in this paper can be compared 
(see e.g. Ref. 0). 

Acknowledgements. The research of TKG and TS is supported in part by the U.S. 
Department of Energy under DE-FG02-91ER40626.A007. The research of SM is supported 
in part by the Grant-in- Aid for Scientific Research, of the Ministry of Education, Science 
and Culture of Japan No. 07640419. 



10 



REFERENCES 

[1] R. Becker-Szendy et al. (1MB Collaboration), Phys. Rev. D 46 (1992) 3720. 

[2] K. S.Hirata et al. (Kamiokande Collaboration), Phys. Letters B 280 (1992) 146 and 
Y.Fukuda et al, Phys. Letters B 335 (1994) 237. 

[3] T.K.Gaisser, Phil. Trans. Roy. Soc. Lond. A 346 (1994) 75. 

[4] E. V.Bugaev and V.A.Naumov, Phys. Lett. B 232 (1989) 391. 

[5] G.Barr, T.K.Gaisser and T. Stanev, Phys. Rev. D 39 (1989) 3532. 

[6] M.Honda, K. Kasahara, K.Hidaka and S. Midorikawa, Phys. Lett. B 248 (1990) 193. 

[7] H. Lee and Y. S. Koh, Nuovo Cimento B 105 (1990) 883. 

[8] W.A.Mann, T.Kafka and W. Leeson, Phys. Lett. B 291 (1992) 200. 

[9] G.L.Fogh, E.Lisi and D. Montanino, Phys. Rev. D49 (1994) 3626. 

[10] O.G.Ryazhskaya, Nuovo Cimento 18 C (1995) 77; Pis'ma Zh. Eksp. Teor. Fiz. 61 (1995) 
229 [JETP Lett. 61 (1995) 237]. 

[11] H.Lee and S. A. Bludman, Phys. Rev. D 37 (1988) 122. 

[12] T. K. Gaisser and T. Stanev, in Proceedings of the 24th International Cosmic Ray Con- 
ference, Rome, Italy, 28 August - 8 September, 1995, Vol. 1, p. 694. 

[13] D.Cooke, Phys. Rev. Lett. 51 (1983) 320. 

[14] T. K. Gaisser, R. J. Protheroe and T. Stanev, in Proceedings of the 18th International 
Cosmic Ray Conference, Bangalore, India, 1983, edited by N. Durgaprasad et al. (Tata 
Institute of Fundamental Research, Colaba, Bombay, 1983), Vol.5, p. 174. (Fig. 1 in this 
reference is incorrectly plotted.) 

[15] L. I. Dorman, V. S. Smirnov, and M. I. Tyasto, Cosmic Rays in the Magnetic Field of the 
Earth, Moscow "Nauka", 1973 (in Russian). 

11 



[16] E. V. Bugaev and V. A. Naumov, Institute for Nuclear Research (USSR Academy of 
Sciences) Report 11-0537 (Moscow, 1987), unpublished. 

[17] V. A. Naumov, in Proceedings of the International Workshop on v^/ve Problem in Atmo- 
spheric Neutrinos, Gran Sasso, Italy, 5-6 March, 1993, edited by V. S. Berezinsky and G. 
Fiorentini (LNGS, L'Aquila, Italy, 1993), p. 25. 

[18] A. N. Kalinovsky, N. V. Mokhov and Yu. P. Nikitin, Transport of High-Energy Particles 
Through Matter, Moscow, "Energoatomizdat" (1985); 

L. R. Kimel' and N. V. Mokhov, University Bulletin, Physics (Tomsk University Press) 10 
(1974) 17; and in Problems of Dosimetry and Radiation Protection, Moscow, "Atomizdat", 
14 (1975) 41; 

A. Ya. Serov and B. S. Sitchev, Transactions of Moscow Radio Engineering Institute 14, 
173 (1973). 

[19] V. A. Naumov, in Investigations in Geomagnetism Aeronomy and Solar Physics, Moscow, 
"Nauka", 69 (Theoretical Physics) (1984) 82. 

[20] E.V. Bugaev, and V. A. Naumov, Yad. Fiz. 45 (1987) 1380 [Sov. J. Nucl. Phys. 45 (1987) 
857]. 

[21] P. Lipari and T. Stanev, in Proceedings of the 24th International Cosmic Ray Conference 
(Ref. [0), Vol.1, p. 516. 

[22] B. Nilsson-Almqvist and E. Stenlund, Comput. Phys. Commun. 43 (1987) 387. 

[23] T. Sjostrand and M. Bengtsson, Comput. Phys. Commun. 43 (1987) 367. 

[24] K. Hanssgen and J. Ranft, Comput. Phys. Commun. 39 (1986) 37. 

[25] L. V. Volkova, in Cosmic Gamma Rays, Neutrinos and Related Astrophysics, edited by 
M.M.Shapiro and J.P.Wefel, NATO ASI 270 (1989) 139. 

[26] V. A. Naumov and A. E. Rjechitsky, Irkutsk State University, Internal Note 3-91 (1991) 



12 



(unpublished) . 

[27] J.Engel, T. K. Gaisser, P. Lipari and T. Stanev, Phys. Rev. D 46 (1992) 5013. 

[28] J. V. Allaby et al, CERN Yellow Report 70-12 (unpublished). 

[29] T.Eichten et al, Nucl. Phys. B 44 (1972) 333. 

[30] T. K. Gaisser, T. Stanev and G. Barr, Phys. Rev. D 38 (1988) 85. 

[31] M.Antinucci et al., Lett, al Nuovo Cimento 6 (1973) 121. 

[32] J.R. Sanford and C.L.Wang, Brookhaven Report BNL 11299/11479 (1967) (unpub- 
lished) . 

[33] S. J.Barish et al, Phys. Rev. D 16 (1977) 3103. 
[34] H. J.Miick et al, Phys. Lett. B 39 (1972) 303. 
[35] D. H. Perkins, Astroparticle Physics 2 (1994) 249. 

[36] G. Basini et al (MASS Collaboration), in Proceedings of the 24th International Cos- 
mic Ray Conference (Ref. |T^) Vol.1, p. 585; J. F. Krizmanic et al (IMAX Collabora- 
tion), ibid, Vol.1, p. 593; E.Schneider et al (HEAT Collaboration), ibid, Vol.1, p. 690; 
R.Bellotti et al (MASS Collaboration), Phys. Rev. D 53 (1996) 35. 

[37] M.Honda, T. Kajita, K. Kasahara and S. Midorikawa, Phys. Rev. D 52 (1995) 4985. 



13 



TABLES 



TABLE L Comparison of calculated neutrinos fluxes at Kamioka 
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TABLE IL Comparison of differential spectra (solar minimum) used by BGS and BN. 



(m ^sr ^s ^GeV Energy is total energy per nucleon.) 
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TABLE III. Comparison of differential spectra (solar average) used by BGS and HKHM. 
(m~^sr~"^s"~^GeV~-^; Energy is total energy per nucleon.) 
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FIGURES 

FIG. 1. Z-factors for charged pions in proton-air collisions. Solid: BGS dotted: HKHM 
dashed: BN @. 

FIG. 2. Ratio of Z-factors for vr^ and vr^. (Same codes as Fig. |^). 

FIG. 3. Average number of charged pions produced by a nucleon with total energy Eq incident 
on air (multiplied by Eq^'"^). Results are shown for three different bins of pion momentum for BGS 
1^] (solid lines) and for BN (dashed lines). 

FIG. 4. Average multiplicity of charged pions per interaction of proton {Ep = total energy) in 



air (except for line which is for proton-proton collisions from Ref. |31]) 



FIG. 5. Distributions of fractional momentum (dn / dlnx) of charged pions produced in inter- 
actions of ~ 20 GeV/c momentum protons with light nuclei. Models are shown for target = air, 



data [28,2^ for target = Be. 7r+ and vr are shown separately. 
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